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ABSTRACT 


Efforts  to  improve  timing  resolution  from  various  scintillator- 

photomultiplier  combinations  are  described  and  the  results  of  this 

work  are  applied  to  the  direct  electronic  measurement  of  nuclear  life- 

45  45 

times  for  levels  of  Sc  and  Ti.  A  lifetime  of  565  +  20  psec  has 

45 

been  found  in  Sc  from  a  level  between  2.8  and  3.2  MeV  excitation 
energy  excited  by  the  ^Sc(p,p')^Sc*  reaction.  In  addition,  upper 

limits  of  180  psec  have  been  set  on  the  mean  lifetimes  of  levels  of 

45  45 

Ti  up  to  1.52  MeV  excitation  energy  strongly  excited  by  the  Sc 

45  * 

(p,n)  Ti  reaction. 

A  review  of  various  lifetime  measurement  techniques  and  improve¬ 
ments  to  the  technique  are  also  included. 
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I .  INTRODUCTION 


This  thesis  will  describe  work  done  on  improving  the  technique  of 
measuring  nuclear  lifetimes  in  the  nanosecond  and  subnanosecond  range  by 
a  direct  electronic  method.  The  method,  a  form  of  delayed  coincidence, 
is  described  in  detail  in  section  III,  but  will  be  briefly  explained 
here  to  provide  a  degree  of  familiarity  with  the  material. 

A  nuclear  state  is  excited  by  regularly  spaced  beam  pulses  from  a 
Van  de  Graaff  accelerator.  The  radioactive  decay  spectrum  is  obtained  by 
measuring  the  time  interval  between  the  beam  pulses  and  the  output  pulses 
from  a  scintillation  detector.  Then,  from  the  radioactive  decay  law 

f  -  -XN  [1] 

where  N  is  the  number  of  still  excited  nuclei  surviving  after  the  time  t 
and  X  is  the  decay  constant,  the  mean  life  of  the  level  is  found 

T  =  y-  and  N  ■  N.e  [2. a] 

m  A  1 

where  N_^  is  the  number  of  excited  nuclei  at  time  t  =  0. 

When  the  number  of  gamma  rays  that  enter  a  detector  in  a  given 
time  period  At  is  set  equal  to  Nq  and  the  number  that  enters  the  same 
detector  under  the  same  conditions  (same  At,  same  experimental  setup)  at 
a  time  t  later  is  set  equal  to  N,  in  the  limit  that  At  *>  0,  the  mean 
life  of  the  gamma  emitting  level  is  given  by 

_ t _ 

Tm  ZnN  -  In  N 
o 

Tm  is  of  interest  to  nuclear  physics  as  its  theoretical  value  is 
extremely  sensitive  to  the  nuclear  model  used  to  calculate  it.  Thus 


T  alii  a  i  s >d3  3£»  ori  ej  A  mb 


/ 
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comparisons  of  theoretical  and  experimental  values  will  provide  sensi¬ 
tive  tests  of  nuclear  models.  This  will  be  shown  in  section  II. 

This  thesis  will  be  divided  into  seven  sections.  The  next  sec¬ 
tion  presents  a  review  of  the  theory  of  nuclear  decay  as  applied  to 
gamma  radiation  and  the  third  section  gives  a  review  and  a  comparison  of 
the  various  methods  of  measuring  nuclear  lifetimes.  The  fourth  section 
describes  the  experiments  carried  out  as  part  of  this  project,  both  on 
the  bench  and  using  the  accelerator,  while  the  fifth  section  contains 
conclusions  and  discussion.  Appendices  are  given  in  section  six.  The 
final  section  will  consist  of  a  glossary.  A  list  of  the  references 


follows  the  final  section. 
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II.  THEORY 


An  excited  nucleus  can  de-excite  by  particle  emission,  8  decay, 
gamma  emission,  internal  conversion  or  fission.  It  is  the  third  method 
which  concerns  us  here  and  this  section  will  present  a  review  of  the 
theory  of  electromagnetic  radiation  with  particular  emphasis  on  the 
transition  probability  which  is  inversely  proportional  to  the  nuclear 
lifetimes.  This  is  done  to  provide  the  general  background  material 
relevant  to  this  topic  as  well  as  to  specifically  relate  the  lifetime 
information  to  nuclear  structure  theory. 

Electromagnetic  radiation  is  classed  by  multipole  order  L  where 
L  is  the  angular  momentum  in  units  of  -ft  carried  off  by  each  emitted 
quantum.  Each  multipole  order  has  two  classes  of  radiation,  the  electric 
2^  pole  (EL)  and  the  magnetic  2^  pole  (ML) .  While  these  differ  with  re¬ 
spect  to  parity  and  classically  they  refer  to  radiation  emitted  by  elec¬ 
tric  or  magnetic  2L  poles,  the  physical  difference  between  electric  and 
magnetic  multipole  fields  can  be  deduced  from  the  following  equation 

r*EM  =  r*HE  =  0  [3] 

which  is  derived  in  Appendix  I.1)*  This  equation  states  that  the  elec¬ 
tric  multipole  field  has  no  non— vanishing  radial  components  of  the  mag¬ 
netic  field  H  (the  electric  multipole  field  does  have  non-vanishing 

*In  equation  3  the  superscript  E  and  M  refer  to  the  electric  or  magnetic 
nature  of  the  radiation  while  the  subscript  M  is  the  Z  component  of  the 
angular  momentum,  L. 
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radial  components  of  the  electric  field  E) ,  while  the  magnetic  multipole 
field  has  non-vanishing  radial  components  of  H.  but  not  of  E. 


The  conservation  of  angular  momentum  and  parity  of  the  nuclear 
system  plus  y-rays  requires 


J.  -  J J  <  L  <  J.  +  J 
if1  1 


[4] 


and 


Ait  =  tt./tt-  =  (-1) 
if 


for  EL  radiation 


[5. a] 


=  (-1)L  1  for  ML  radiation  [5,b] 

where  J_^  and  are  the  initial  and  final  angular  momenta  of  the  nuclear 
states  (parities  tt^  ,  tt^),  and  L  is  the  angular  momentum  carried  off  by 
the  emitted  quantum.  Att  =  +1  denotes  no  parity  change  while  Air  =  -1 
denotes  a  change  of  parity.  In  addition,  no  radiation  of  multipole 
order  0  can  be  emitted.  This  is  a  consequence  of  the  transverse  nature 
of  electromagnetic  radiation  and  can  be  seen  from  equation  I.2l  in  Ap¬ 
pendix  I. 


In  a  study  of  the  gamma  radiation  transition  probability  the 
single  particle  model  is  first  considered.  In  particular,  the  transi¬ 
tion  probability  for  a  single  non-relativistic  proton  moving  in  a 
potential  due  to  all  the  other  particles  in  the  nucleus  is  evaluated. 
This  is  not  a  correct  picture  of  the  nucleus,  but  as  transition  proba¬ 
bilities  of  more  realistic  models  can  be  expressed  in  terms  of  the 
results  of  this  simpler  view,  we  proceed. 

The  quantum  mechanical  transition  probability  T  f  for  any  pro¬ 
cess  is  given  by 
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2tt 

ft 


<f |tf' | i> | 


2  dN 
dE 


[6] 


where  H'  is  the  perturbing  interaction  and  |i>  is  the  first  order 

matrix  element  of  the  interaction  between  the  wave-function  of  the  ini¬ 
tial  state  i  and  the  final  state  f.  denotes  the  density  of  final 

states.  This  equation  is  developed  in  Appendix  II.  In  a  non-relativ- 
istic  treatment  the  total  Hamiltonian  of  the  proton  and  the  electromag¬ 
netic  field  is 


H  = 


(P_  “  (e/c)A) 
2m 


+  V (r)  -  y  a/H  +  £  n  hu) 


p  2mc 


X  X 


[7] 


where  m  is  the  proton  mass,  y^(efi/2mc)  its  magnetic  moment,  H  is  the 
magnetic  field,  o_  the  Pauli  spin  vector,  and  P_  is  the  momentum  operator 
-ift  grad.  In  the  lowest  order  of  interaction  (single  quantum  emission) , 
this  becomes 


H' (A)  =  - 


P • A  +  A*P 
e - - 

me  2 


eh 

2mc 


cr-H 


[8] 


The  Hamiltonian  is  taken  between  an  initial  state  consisting  of  a  pro¬ 
ton  in  state  i  and  no  quanta  and  a  final  state  f  with  one  quantum  char¬ 
acterized  by  quantum  numbers  X.  From  equations  1.4  and  I. 11  of  Appendix 


I 


where 


<f  ,1A  | H'  (A)  I  i  ,0> 


*1/2  * 

<fl*'<Vli> 


[9] 


[10] 


Note  that  the  matrix  element  on  the  right  hand  side  of  equation  9  is 

independent  of  q  ,q,  (the  creation  and  annihilation  operators  of  Ap- 

A  A 


pendix  I),  and  time. 


. 
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The  transition  probability  term  for  the  density  of  states  in  a 
spherical  volume  of  large  radius  Ro  is 


dN  _  Ro 
dE  Trftc 


[11] 


Thus 


[12] 


Applying  equations  1.22  to  equation  8 


=  “  ^-I'(curl  k  «„,) 


*  CT  y_efi 


Lp  „ 

-  “2  *£•  — ^LM 

2mc 


[13. a] 


M*  C  ieP-(Lu  )*  C  tii  eft 

*'  ^  •  -S£-  &•  aw 


[13. b] 


At  this  point  the  approximation  that  the  nuclear  radius  a  is  small  com¬ 
pared  to  the  wavelength  of  the  radiation, 


«  1 
c 


[14] 


is  introduced,  allowing 


„  =  (ur/c)  y  aN 

ULM  (2L+1) ! !  ^LM  ^ 6  *  ^ 


[15] 


(from  equation  11.15  and 


(kAr) 

^L(kAr)  k^r+0  (2L+1)  !  !  ^ 


[16] 


Stech  )  has  proven  the  following  two  relations  which  are  used  to  further 
simplify  the  expressions  for  the  transition  probabilities: 


(curl  Lu^)*  =  ^[£AuLh  x  grad  +  2)8rad  u^]*  [17. a] 


-  i  (L+l)  [grad  u^]* 


[17. b] 


and 
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<f  IJP-  (grad  uLM)  |i>  =  - 


2  2 
.  iP  P  |  .  v 

*  <f  2m  “LM  °LM  2m‘1> 


=  -  inuo<f  |  u^  |  i> 


[18] 


Using  these  and  introducing  the  quantities  (a  5  E  or  M) 


=  CL  "(2L+jyTT  (|)L<f|<„|i>  [19] 


L+l  ,U)vL 


where 


—  er^y 
LM 


*  in  ehiD  ,  T  * 

LM  -  -f-T"  (L+1)  2.x£-[grad(rLy  )] 
2mc 


[20. a] 


^LM  =  mc^  L+l  — '  tgrad(rLJLM)]*  +  £  I*  a-  [grad<rL7  )]*  [20.b] 


Using  the  expression  for  C  from  equation  1.25 

Li 


T(aLM) - 8tt  (L+l)  1  ^  2L+1 1  <f  ,.>|2 

L[(2L+l)!!r  11  C  LM 


i-*f 


[21] 


Before  proceeding  to  derive  the  single  particle  transition 
rates  (in  terms  of  which  the  more  realistic  rates  are  expressed)  several 
observations  may  be  made.  While  the  selection  rules  given  by  equations 
4  and  5  still  hold,  equations  20  imply  an  additional  criterion  on  the 
orbital  angular  momenta  (Z^,Zf)  if  the  particle  is  moving  in  a  central 
velocity-independent  potential, 

\l.  -  l  |  <  L  <  l±  +  lf  [22] 

Equations  20  and  21  can  also  be  used  to  show  that  if  (wa/c)<<l,  then 
transitions  for  a  single  proton  will  tend  to  proceed  by  the  lowest  mul¬ 
tiple  order  as  follows.  By  taking  the  orders  of  magnitude  for  the  A^M, 

*  o(ea^) 


[23. a] 


brrs 


fiji  .act  x  ,.fu  ..  IfiJldio 
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■  °t(^7HeaL>]  •  [23. b] 

orders  of  magnitude  for  the  transition  probabilities  can  be  expressed  as 
follows 


2 

=  °tu)^)(a>a/c)2L]  t24-ai 

2 

t(ML)  _  0[o)(~^)  (*/mca)2(a)a/c)2L]  [24. b] 


The  ratios  of  the  two  lowest  terms  are: 


1)  |ji  -  jf|  -  L(i  1),  Air  =»  (-1)L 


2) 


(M(L+1)) 

- —  =  o[(ho)/mc  )  ] 


|j±  -  jf|  =  L(2  1),  Air  =  (-1)L_1 


[25. a] 


(E(L+1)) 

-^THE5 - °t("a2/*)2] 


[  25 .  b  ] 


In  general  the  ratios  of  equations  25  will  be  much  less  than  unity  and 
the  radiative  transitions  will  in  fact  proceed  by  the  lowest  multipole 
order.  In  actual  practice,  many  transitions  involve  the  two  lowest 
orders  and  more  complex  models  do,  in  fact,  account  for  this.  It  is  in 
deducing  the  single  particle  transition  probabilities  that  the  sensi¬ 
tivity  of  the  study  of  nuclear  lifetimes  to  nuclear  structure  becomes 
evident.  The  dependence  of  the  transition  probability  on  the  form  of 
the  initial  and  final  state  wavefunctions  makes  the  study  of  nuclear 
lifetimes  an  excellent  check  on  the  validity  of  these  model-dependent 
wavefunctions.  This  is  seen  as  follows.  If  the  proton  is  assumed  to  be 


ovi  sril  sv  o'  nl  an  11  an  nr  he  ,  >J  t:  .  enl  >c  a  .sabio 


moving  in  a  central  potential  each  wavefunction  and  can  be  written 
as  a  product  of  a  radial  function  and  an  angular  function.  The  radial 
functions  R^  and  R ^  depend  upon  the  details  of  the  potential,  while  the 

mi  mf 

two  component  functions  07  .  and  07  .  contain  angular  and  spin  vari- 

i » J  i  f  f 

ables  only.  This  allows  the  electric  multipole  matrix  element  to  be 
rewritten 


<^1°  ■ 


roo 


R  r^R . r2dr 
f  1 


m. 


m 


V  .  y* e,1  .  an 

4.  lvh  LM  h^i 


[26. a] 


where  the  contribution  from  the  spin  term  in  equation  20. a  has  been 
neglected  as  insignificant.  However,  the  magnetic  multipole  matrix  ele¬ 
ment  is  much  more  complex,  as  the  spin  term  is  actually  more  important 
than  the  orbital  term  and  so  a  simplifying  special  case  is  considered, 

| j -  jf|  =  L  and  \l  -  l  \  =  L  -  1.  With  these  criteria,  equation 
20. b  gives 


=  2^  %  -  LfT)<f|H.-8rad<rLyM>li>  t26-b] 


The  transition  probability  is  summed  over  the  m  substates  of  the 
final  state  and  is  averaged  over  the  m  substates  of  the  initial  state 
to  give  a  total  transition  probability 


-  (2ji  +  i)'l  l  It 


mi  mf  ^ 


,(oLM) 

i+f 


[27] 


This  is  evaluated  with  the  help  of  a  statistical  factor 


S(j.,L,jf)  =  4ir(2j  +  1)  1  I  I  II 
11  1  m.  m_  M 


i  f 


mf  *  ,2 

0/  .  7tm07  .  dft  [28. a] 

4*  VJf  Zi’di 


-  (2j  +  l)[C(j.,jf,L, 1/2, -1/2,0)] 


[28. b] 


.  (L  5  3  II  X  ".0  ow3 


oM  to  as^£3i  i  m  sri  i  isvo  bommua  i  ^  ifioo'rq  nol  ansii  orf? 
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where  C  is  a  Clebsh-Gordan  Coefficient.  The  values  of  S(ji>L,jf)  are 
listed  in  tables.'*')  Using  this  statistical  factor,  the  transition  prob¬ 
abilities  become 


(EL) 

i+f 


t(ML) 

i+f 


2 (L+l)  ,e\  ,  ,  .2L, 

- — 7  o)(~)  (oaa/c)  ( 

hi2  ftc 


L[(2L+1)  !  !  ] 


2  (L+l)  ,e*  .  v  2L 

-  co (— )  (u) a/c) 


L[ (2L+1) ! ! ] 


til2 


rco 


y  O 
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By  noting  that  =  —  ,  sensitive  tests  of  the  form  of  the  R's  may 

m 

be  made  providing  clues  to  nuclear  structure.  If  the  radial  functions 
are  assumed  to  be  constant  inside  the  nuclear  volume  (r  <  a)  and  vanish 
outside  (r  >  a) , 

f00 

(r/a)IjR^r2dr  =  3(L  +  3)  ^  .  [30] 

o 

With  this  assumption,  the  single  particle  transition  probabilities  for 
a  single  proton  become 
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where  ftu>  is  a  MeV,  and  a  is  in  fermis.  Table  I  gives  the  estimates  for 
various  multipolarities. 
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TABLE  I 


Multipolarity 

T 

sp 

,  -1 
in  sec 

El 
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a2/3 

Ey 

Ml 

13 

2.9  x  10 

Ey 

E2 

7.4  x  107 

a4/3 

Ey 

M2 

8.4  x  107 

a2/3 

Ey 

E3 

3.4  x  101 

A2 

Ey 

M3 

8.7  x  101 

A4'3 
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E4 

1.1  x  10"5 

a8/3 

Ey 

M4 

4.8  x  10-5 
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For  single  neutron  transitions  the  electric  transition  probabilities  are 
smaller  than  the  proton  transition  probabilities  by  a  factor 

Z/AL  [32. a] 

and  the  magnetic  transition  probabilities  are  slightly  smaller  than  the 
corresponding  probabilities  for  electric  transitions,  the  ratio  being 


[32. b] 


Mention  should  be  made  of  many  particle  configurations  for  which 
a  reduced  transition  probability  is  defined 
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This  gives  for  the  transition  probability 
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III.  REVIEW  OF  METHODS  OF  LIFETIME  MEASUREMENT 

The  purpose  of  this  section  is  to  review  the  methods  used  in  the 

study  of  nuclear  lifetimes.  Such  a  review  has  been  carried  out  by 

3 

Schwarzschild  and  Warburton  )  and  the  techniques  dealt  with  in  their 
paper  are  treated  below  along  with  a  brief  mention  of  Resonance  Absorp¬ 
tion  and  Fluorescence.  A  division  is  made  between  methods  employing  a 
form  of  Doppler  shift,  and  those  which  do  not.  The  former  class  includes 
the  Doppler  Shift  Attenuation  Method  (DSAM) ,  Resonance  Methods  and  the 
Recoil-Distance  method.  Delayed  coincidence  methods  including  the  Direct 
Electronic  Timing  method  make  up  the  latter  category. 

The  first  Doppler  shift  technique  is  the  Doppler  Shift  Attenua- 

j\  *  £.  . 

tion  Method  (DSAM).  This  method  has  been  discussed  in  detail  by  &*A. 
Litherland  et  al.S  and  E.K.  Warburton  et  al.5).  It  involves  exciting 
a  target  nucleus  which  particle  decays  to  an  excited  state  that  sub¬ 
sequently  gamma  decays.  After  the  particle  decay,  the  excited  residual 
nucleus  recoils  in  accordance  with  momentum  conservation  and  is  slowed 
down  at  a  known  rate  in  the  target  backing  material.  The  gamma  rays 
emitted  are  Doppler  shifted  to  a  degree  dependent  upon  the  velocity  of 
the  recoiling  nucleus  at  the  time  of  gamma  emission.  By  measuring  this 
shift  and  by  knowing  the  recoil  velocity  as  a  function  of  time,  the 
time  of  emission  and  hence  the  lifetime  can  be  determined. 

-14 

The  method  is  effective  for  lifetimes  in  the  range  10  to 
10-11  secS  for  solid  targets  and  the  range  is  extended  to  10  10  sec6) 


14 

for  gas  targets.  While  this  implies  that  the  choice  of  target  is  severe¬ 
ly  limited  for  lifetimes  longer  than  10  ^  sec,  the  method  has  the  advan¬ 
tage  of  being  relatively  quick  and  easy  to  carry  out.  However,  the 
precise  knowledge  required  of  the  recoil  ion  velocity  as  a  function  of 
time  limits  the  accuracy  of  the  DSAM  to  the  order  of  15  -  20%^).  Other 
effects  such  as  Doppler  broadening  and  inhomogenieties  in  the  target 
material  and  backings  also  give  rise  to  uncertainties.  But  above  all, 
the  method  while  effective,  is  an  indirect  one  and  it  is  with  this  point 
in  mind  that  attempts  are  being  made  to  extend  the  more  direct  methods 
into  this  lifetime  region. 

The  second  Doppler  shift  method  to  be  studied  is  really  two 
variations  of  the  same  theme,  Resonance  Absorption  and  Resonance  Fluo¬ 
rescence.  In  atomic  physics  resonance  absorption  is  a  well  known  phe¬ 
nomenon.  The  absorption-line  spectrum  where  certain  frequencies  of  an 
impinging  light  beam  excite  the  orbital  electrons  of  a  gaseous  material 
and  are  thereby  lost  to  the  beam  is  an  old  analytical  tool  often  used 
to  identify  a  material  in  question.  A  related  phenomenon  is  fluores¬ 
cence  or  resonance  scattering  in  which  the  impinging  light  causes  res- 

g 

onant  stimulated  emission  of  radiation.  Although  Kuhn  )  predicted  as 
early  as  1929  that  these  resonant  features  would  occur  in  nuclei,  they 

9 

have  been  observed  only  within  the  past  two  decades  initially  by  Moon  ) . 
This  is  because  the  recoil  of  the  emitting  nucleus  takes  energy  and 
Doppler  shifts  the  frequency  of  the  emitted  radiation  into  a  frequency 
which  cannot  be  reabsorbed  by  a  second  similar  medium. 


c 


[35] 
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where  /  is  the  frequency  of  the  radiation  as  seen  by  the  second  nucleus 

which  would  absorb  the  radiation  if  the  shift  were  not  taking  place,  and 

is  tke 

f  *w4th  natural  emission  frequency. 

These  phenomena  can  be  used  to  measure  nuclear  lifetimes  in  the 
following  way.  The  spectral  lines  have  a  natural  width  in  their  frequen¬ 
cies  600  which  is  proportional  to  the  mean  life  t  of  a  nuclear  state  ac¬ 
cording  to  the  relation 

rx  -  *  [36. a] 

and 

T  -  A6o>  [36. b] 

where  T  is  the  level  width  in  energy  units,  and  “ft  is  Planck's  constant 
divided  by  2tt.  Thus  if  the  level  width  can  be  found,  the  lifetimes  of 
the  state  in  question  will  be  known.  Table  II  lists  approximate  values 
of  T  with  corresponding  values  of  x. 


TABLE  II 


r(ev) 

x(sec) 

r(eV) 

x (sec) 

6  x  lcf 10 

10-6 

6  X  10’1 

10" 15 

6  X  10"7 

io“9 

6  X  102 

i-* 

o 

1 

t-> 

00 

6  x  10  ^ 

i<f12 

6  X  105 

rH 

CM 

1 

O 

i— 1 

6co  can  be  found  if  two  pieces  of  the  same  material  are  set  up  so 
that  one  piece  which  is  excited  and  emits  gamma  rays  can  be  made  to  move 
towards  a  non-excited  piece  which  acts  as  a  resonant  scatterer  or  ab¬ 
sorber.  This  motion  of  the  source  medium  can  be  made  fast  enough  mechan- 
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ically  so  that  it  creates  a  Doppler  shift  of  the  emitted  gamma  ray  which 
will  be  exactly  equal  and  opposite  to  the  recoil  Doppler  shift,  thereby 
allowing  the  'target'  piece  to  act  as  a  resonant  scatterer  or  absorber. 
By  measuring  the  intensity  of  the  scattering  or  absorbing  as  a  function 
of  the  source  medium  velocity  (hence  as  a  function  of  energy  due  to 
Doppler  shifting),  the  degree  of  the  overlap  of  the  shifted  gamma  ray's 
natural  energy  spread  and  the  energy  level  in  question  (the  scattering 
or  absorption  cross  section)  gives  the  width  of  the  level  in  terms  of 
energy  units.  In  practice  the  method  is  limited  in  several  ways.  At 
the  outset,  the  Doppler  shift  energies  must  be  resolveable  implying  that 
this  method  is  applicable  only  to  nuclear  lifetimes  in  the  10  ^  sec 
region  and  shorter  (see  table  II).  Also,  the  isotope  in  question  must 
not  only  have  a  substantial  natural  abundance  to  be  useful  as  a  scat¬ 
terer,  but  it  must  also  be  available  in  an  excited  state  following  the 
decay  of  its  radioactive  parent.  One  of  the  more  formidable  drawbacks 
is  that  the  velocity  required  for  the  correct  Doppler  shift  be  readily 
attainable  by  mechanical  means  (up  to  10^  cm/sec).  These  considerable 
problems  are  in  most  cases  insurmountable  leaving  the  technique,  at 
best,  a  demonstration  of  the  correctness  of  the  theory  of  resonance 
fluorescence. 

In  1958,  Mossbauer10)  discovered  another  method  of  applying  res¬ 
onance  effects  to  lifetimes  studies,  in  this  case  recoilless  resonance 
absorption  rather  than  the  resonance  fluorescence  specifically  studied 
by  Moon.  Mossbauer  reasoned  that  if  the  emitting  nucleus  had  a  very 
large  mass  it  would  receive  the  recoil  momentum  but  almost  no  energy 
from  a  gamma  transition,  thus  almost  completely  eliminating  the  gamma 
ray  Doppler  shift.  In  fact,  an  infinitely  heavy  nucleus  does  not  exist. 
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but  a  nucleus  bound  in  a  rigid  crystal  will  act  as  an  infinite  mass  and 
recoilless  transitions  can  thus  take  place.  It  should  be  noted  that  the 
term  ' recoilless ’  refers  to  energy  only  and  not  to  the  momentum  trans¬ 
ferred  to  the  emitting  nucleus.  Under  these  conditions,  small  relative 
velocities  between  the  source  and  absorber  bodies  may  be  used  to  measure 
the  absorption  cross  section  giving  the  level  width  and  hence  the  life¬ 
time  involved.  While  this  eliminated  the  high  velocity  problem  of  Moon’s 
technique,  the  other  drawbacks  mentioned  above  still  apply. 

11-14 

In  addition,  a  detailed  account  of  the  theory  of  this  method  ) 

introduces  an  important  factor,  the  Debye-Waller  factor,  into  the  final 

15  16 

result  and  this  can  be  measured  to  within  not  better  than  20%  )  ). 

The  cryogenic  temperatures  required  for  the  experiment  and  accurate  low 
velocity  drives  are  two  more  features  which  make  this  method  difficult 
to  use  for  lifetime  studies.  As  a  result,  extensive  work  has  not  been 
carried  out  with  this  method. 

The  final  Doppler  shift  technique  to  be  reviewed  is  the  Recoil- 

Distance  method.  It  is  intended  to  be  used  in  the  region  between  the 

-9 

DSAM  region  and  the  direct  electronic  timing  region,  i.e.  from  10  to 
5  x  10  seconds.  The  technique  is  described  in  detail  by  Alexander 
and  Allen  )  but  is  briefly  presented  as  follows.  A  heavy  ion  beam  is 
used  to  strike  a  thin  target  back-mounted  on  a  thin  supporting  foil 
causing  the  residual  nuclei  to  recoil  freely  into  a  vacuum  with  velocity 
V.  Since  the  characteristic  slowing  down  time  for  heavy  ions  in  a  metal 
is  typically  3  x  10  sec,  a  metal  plate  on  a  moveable  plunger  placed 
at  a  distance  s  =  vt  from  the  target  will  stop  the  recoiling  ions  in  a 
time  which  is  very  short  compared  to  their  nuclear  lifetimes.  The  in- 
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tensity  of  the  gamma  rays  from  nuclei  slowed  down  in  the  plate  will  thus 

give  the  number  of  nuclei  surviving  for  a  time  t.  The  gamma  rays  from 

nuclei  which  are  stopped  (i.e.  ions  that  have  reached  the  metal  plate 

with  lifetimes  >  t)  are  separated  from  those  which  decay  in  flight  by 

using  their  Doppler  shift  separation  in  energy.  By  varying  the  distance 

s,  a  decay  spectrum  is  built  up  as  a  function  of  t.  While  this  method 

18 

can  use  a  variety  of  solid  (and  gas  ))  targets  in  a  fairly  wide  time 
range,  not  only  does  it  suffer  from  being  an  indirect  method  but  only 
heavy  ions  can  be  used  as  projectiles  and  Doppler  broadening  of  the 
shifted  line  reduces  accuracies  to  about  10%^). 

The  other  basic  method  of  measuring  nuclear  lifetimes  is  by  the 

delayed  coincidence  technique  and  its  variations,  of  which  our  method 

is  one.  In  its  simplest  form  a  radioactive  parent  nucleus  with  a  daugh- 

-7  -10 

ter  possessing  a  mean  gamma  decay  lifetime  in  the  range  10  to  10 

sec  is  considered.  Two  detectors  are  used,  one  starting  a  clock  at  the 

time  of  the  creation  of  the  daughter  nucleus  by  recording  an  emitted 

particle  from  decay  of  the  parent  nucleus  and  the  second  stopping  the 

clock  by  recording  the  gamma  ray  from  the  decay  of  the  daughter  nucleus. 

By  taking  many  of  these  measurements  an  exponential  decay  spectrum  is 

built  up  of  gamma  rays  v.s.  time,  and  the  mean  life  t  is  found.  While 

this  technique  is  standard  and  has  been  used  a  great  number  of  times 
19 

(cf.  ))  it  is  limited  in  the  choice  of  nuclei  to  be  studied  by  the 

parent-daughter  requirement.  It  is  in  overcoming  this  that  the  direct 

20 

electronic  method  was  developed  ) .  A  pulsed  beam  is  made  to  strike  a 

1  * 

target  composed  of  either  the  nuclei  to  be  studied  (M^Cn^  ,m2)M3  reac¬ 
tions)  or  a  nucleus  which  can  be  transmutted  into  the  required  nucleus 
(M1(m2,m3)M*  reactions).  The  time  of  excitation  of  the  final  nucleus  is 
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given  by  having  the  arrival  time  of  a  beam  pulse  on  the  target  start  a 

clock,  while  the  emitted  gamma  ray  is  used  to  stop  the  clock.  After  many 

beam  pulses  an  exponential  decay  spectrum  is  collected  and  this  is 

analysed  to  give  t,  the  lifetime.  Analysis  is  generally  carried  out  by 

two  methods,  least  square  fitting  to  either  the  exponential  decay  curve, 

and  centroid  shift  (see  below) .  These  methods  are  quite  simple  when 

only  one  gamma  ray  with  a  lifetime  greater  than  the  system  resolution 

(cf.  section  IV)  is  present.  Both  methods  require  a  ’prompt*  peak  to 

be  taken,  i.e.  a  spectrum  of  gamma  rays  from  levels  whose  lifetimes  are 

so  short  that  their  exponential  decay  tails  are  completely  masked  by  the 

resolution  of  the  detector  system.  The  centroid  shift  method  due  to 
21 

Bay  )  involves  the  moments  of  the  distribution  of  the  gamma  rays  with 

th 

respect  to  time  where  the  n  moment  is  given  by 


PKAJ 

M^n^[f(x)]  =  |  xnf(x)dx 


[37] 


Bay  considered  only  the  first  moment  to  calculate  the  lifetime 


T  = 


M(1) [N  (t) ] 

_ S _ 


M(1)[v  (t)] 

_ S _ 


[38] 


where  N  (t)  is  the  delay  curve  distribution  and  v  (t)  is  the  prompt  curve 
g  g 

distribution  where  both  are  normalized  such  that 

•OO  f 00 


N  (t)dt  = 


v  (t)dt  =  I  , 


[39] 


I  thus  being  the  area  under  the  curves.  Holland  and  Lynch  )  tried  using 
higher  moments  as  well  but  found  that  less  intense  lifetime  contributions 
from  other  gamma  rays  tended  to  disrupt  the  calculations. 


The  least  squares  method  involves  fitting  either  a  calculated  ex- 
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ponential  curve  to  the  decay  spectrum  to  find  values  of  Nq  and  X  (equa¬ 
tion  2. a)  or  fitting  a  straight  line  to  the  natural  logarithm  of  the 
curve  (equation  2.b).  In  either  case  it  is  usual  for  the  prompt  peak  to 
be  subtracted  from  the  data  so  that  any  slope  present  in  the  logarithm 
plot  is  due  to  the  ’lifetime’  gamma  ray  only. 

The  lifetime  region  of  this  method  can  be  extended  in  the  ’long 
time  length'  direction  almost  indefinitely  and  is  limited  in  the  short 
lifetime  direction  by  the  fact  that  the  decay  tail  must  have  a  slope 
which  is  less  than  that  of  the  system  resolution  peak  if  the  decay  tail 
is  to  be  seen.  Thus  a  great  deal  of  effort  has  gone  into  making  the 
system  resolution  or  prompt  peak  slope  as  narrow  and  as  steep-sided  as 
is  possible  (cf.  section  IV).  The  limit  of  our  work  along  this  line  is 
now  in  the  10  ^  sec  region. 

Obvious  advantages  of  this  method  are  the  variety  of  target  ma¬ 
terials  and  constructions  allowed  and  the  fact  that  this  method,  more 
than  any  other,  is  a  direct  measurement  of  the  nuclear  lifetimes. 
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IV.  EXPERIMENTAL 


A  simplified  diagram  of  the  experimental  setup  is  given  in 

figure  l.a  and  is  explained  as  follows.  The  University  of  Alberta  Van 

de  Graaff  accelerator  is  used  to  provide  a  pulsed  beam  of  protons  at 

1  MIhz  and  0.4  nsec  by  a  Mobley  bunching  system.  This  system  has  been 

22 

described  in  detail  in  previous  publications  ) .  A  pulse  of  the  beam 
first  passes  through  the  capacitive  pickoff  to  produce  a  signal  which 
is  delayed  and  used  to  create  the  stop  signal  for  the  timing  clock. 

The  beam  pulse  then  strikes  and  excites  the  target  which  emits  a  gamma 
ray  after  a  time  characteristic  of  the  target  material. 

In  the  experiment,  gamma  radiation  is  detected  with  scintilla¬ 
tors  composed  of  a  phosphorescent  material  in  a  supporting  medium.  In 
the  scintillator  the  incident  gamma  ray  Compton  scatters  on  an  electron 
and  the  latter  slows  down  by  exciting  the  atoms  of  the  medium  and  the 
phosphorescent  material.  Thus  the  maximum  amount  of  light  emitted  for 
an  incident  gamma  ray  is  proportional  to  the  maximum  Compton  recoil 
energy.  This  emitted  light  is  transformed  into  an  electronic  pulse  by 
a  photomultiplier  tube  which  consists  of  a  sensitive  photocathode 
(which  gives  off  electons  when  struck  by  the  light)  followed  by  electron 
multiplication  stages.  Thus  when  a  gamma  ray  strikes  the  scintillator 
a  pulse  is  created  which  is  used  to  feed  the  fast  start  discriminator 

which  in  turn  produces  a  signal  used  to  start  the  clock.  The  size  of 

lnvers«ly 

the  clock  output  signal  is  proportional  to  the  length  of  the  delay 

't 

between  excitation  and  decay  of  the  target  nuclei  (the  nuclear  lifetime) 
and  this  analogue  signal  is  converted  into  digital  form  which  the 
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computer  can  process  by  the  Analogue  to  Digital  Converter  (ADC).  The 


computer  is  used  to  accumulate  and  display  many  of  the  events  in  a 
decay  spectrum  whose  precise  shape  is  a  function  of  the  nuclear  lifetime 
involved.  This  circuit  is  generally  run  with  a  neutron  elimination 
circuit  (figure  l.b)  to  separate  events  from  the  gamma  rays  from 
spurious  neutron  induced  events.  If  one  and  only  one  gamma  ray  with  a 
lifetime  in  the  region  of  interest  (termed  a  "long  lived"  gamma  ray) 
is  produced  in  the  interaction  of  the  beam  with  the  target,  the  lifetime 
of  the  level  producing  this  gamma  ray  is  given  by  the  natural  logarithm 
of  the  slope  of  accumulated  decay  curve  (cf.  equation  2).  However,  when 
more  than  one  such  gamma  ray  is  present,  or  when  the  identity  of  the 
level  producing  a  given  gamma  ray  is  in  question,  identification  tech¬ 
niques  must  be  introduced.  One  useful  technique  is  to  record  a  Ge(Li) 
spectrum  of  the  gamma  radiation  as  the  mode  of  data  collection  described 
just  above  does  not  differentiate  between  gamma  rays  of  different  energy 
or  origin.  The  Ge(Li)  spectrum  will  allow  such  a  determination  to  be 
accurately  made.  A  further  sophistication  of  this  technique  involves 
gating  the  Ge(Li)  spectrum  with  gamma  rays  which  fall  in  the  "lifetime 
tail"  region  of  the  accumulated  decay  spectrum  so  that  the  Ge(Li) 
spectrum  will  then  show  only  those  gamma  rays  which  are  produced  in 
cascade  (i.e.  simultaneously)  with  the  "lifetime"  gamma  rays.  This  can 
help  identify  gamma  rays  from  long  lived  levels  which  appear  in  the 
"lifetime  tail"  of  the  decay  curve.  However,  this  technique  suffers 
predominantly  from  a  low  count  rate  and  so  was  not  used  in  this  study. 

The  technique  that  was  used  makes  use  of  the  dependence  of  the 
scintillation  detector’s  output  pulse  height  with  gamma  ray  energy.  It 
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i$  called  the  two  dimensional  (TWOD)  collection  and  display  mode  and  gives 
timing  information  along  one  axis  and  pulse  height  information  (gamma 
ray  energy)  along  the  second.  The  number  of  counts  in  a  given  channel 
are  displayed  in  the  "third"  dimension. 

For  each  incoming  signal  from  the  detector,  the  pulse  height 
information  from  the  slow  output  of  the  detector  is  used  to  sort  the 
timing  signals  into  bins,  the  lower  pulse  height  signals  going  into 
lower  bins  and  the  larger  pulse  height  signals  going  into  higher  num¬ 
bered  bins.  Due  to  the  dependence  of  the  pulse  height  on  gamma  ray 
energy,  the  higher  energy  gamma  rays  give  decay  spectra  in  the  upper 
bins  (and  the  lower  bins  too,  due  to  Compton  scattering  in  the  detector) 
while  the  lower  energy  gamma  rays  give  spectra  only  in  the  lower  bins. 
With  a  correct  calibration  of  the  second  analogue  to  digital  converter 
(ADCB)  (cf.  figure  l.b)  which  handles  the  pulse  height  information 
input  to  the  computer,  the  maximum  Compton  transfer  energy  and  hence 
the  full  gamma  energy  can  be  found  for  each  bin. 

Thus,  if  a  decay  curve  is  present  in  the  lower  bins,  by  noting 
the  bin  number  above  which  it  is  no  longer  found  a  rough  indication  of 
the  energy  of  that  gamma  ray  may  be  ascertained.  By  then  referring  to 
the  Ge(Li)  spectrum  mentioned  above,  identification  of  the  gamma  ray 
energy  and  hence  its  level  of  origin  may  be  made. 

Determination  of  the  lifetime  is  made  as  in  the  simple  case  if 
the  bin  in  question  has  only  one  long  lived  contributor.  When  a  given 
bin  has  more  than  one  long  lived  component,  an  "unfolding"  procedure 
must  then  be  used.  Assume,  for  simplicity,  that  th^e  gamma  rays  con¬ 
tribute  to  the  time  spectrum  in  a  given  bin,  two  of  these  being  from 
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long  lived  states  (one  much  longer  (#1)  than  the  other  (#2))  and  one 
(#3)  which  comes  from  a  state  whose  lifetime  is  so  short  that  the 
intrinsic  resolution  of  the  circuitry  masks  its  gamma  decay  curve.  To 
start  the  unfolding  procedure  the  natural  logarithm  of  the  decay  curve 
in  the  bin  in  question  is  plotted.  An  example  is  shown  in  figure  l.c. 

The  straight  line  at  the  left  is  due  to  the  decay  tail  of  the 
longest  lived  state  (#1) .  This  lifetime  may  be  measured  by  applying 
equation  2.b  to  this  line  between  the  points  A  and  B.  Point  B  is 
chosen  so  that  the  contribution  from  the  second  longest  lived  state  (#2) 
has  fallen  to  a  value  which  is  less  than  the  statistical  error  at  that 
point.  The  antilogs  are  taken  of  the  values  on  the  broken  line  between 
B  and  C  and  these  are  then  subtracted  from  the  antilogs  of  the  points 
B  and  E.  When  this  is  done  the  natural  logarithm  is  taken  of  the  new 
curve  between  B  and  E  and  this  is  plotted  again.  This  new  plot  will 
have  a  straight  line  between  F  and  D,  the  latter  chosen  so  that  the 
contribution  from  the  system  resolution  (which  masks  the  shortest  state) 
has  fallen  to  a  value  which  is  less  than  the  statistical  error  at  that 
point.  Applying  equation  2.b  to  this  straight  line  gives  the  lifetime 
of  the  #2  gamma  ray  and  the  whole  process  is  repeated  to  get  the  system 
resolution  slope  which  in  effect  defines  the  shortest  lifetime  which  can 
be  measured.  When  the  unfolding  process  is  completed  (in  this  case 
after  two  subtractions  have  been  made) ,  all  that  is  left  is  a  peak  due 
to  the  system  resolution  called  the  "prompt"  peak.  It  is  the  same 
prompt  peak  that  would  be  displayed  if  there  were  only  very  short  life¬ 
time  gamma  rays  being  collected  in  a  given  bin,  gamma  rays  with  life¬ 
times  so  short  that  the  slopes  of  their  natural  logarithmic  plots  are 
steeper  than  the  slope  due  to  the  system  resolution. 
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From  the  above,  it  can  be  seen  that  if  shorter  lifetimes  are 
to  be  measured  the  system  resolution  or  prompt  peak  should  be  made  as 
narrow  as  possible  to  avoid  masking  the  slope  of  interest  by  the  prompt 
peak.  For  this  reason  a  great  deal  of  effort  went  into  bench  testing  of 
the  detector  system  to  improve  the  overall  resolution.  Before  describ¬ 
ing  in  detail  the  experimental  setup,  the  bench  testing  will  be 
described  to  indicate  the  reasons  for  the  choice  of  the  components 
actually  used. 

Tests  were  made  of  various  photomultiplier  tubes  and  scintilla¬ 
tors  to  try  and  achieve  the  best  possible  timing  characteristics.  At 
the  time  work  began  on  this  thesis  study,  the  best  full  width  at  half 

maximum  (FWHM)  achieved  was  ~210  psec  using  gamma  radiation  from  the 
60 

decay  of  Co  to  provide  coincident  input  to  two  detectors.  The 

circuit  of  figure  2,  with  1"  diameter  cylindrical  Naton  scintillators 

mounted  on  RCA  8575  photomultiplier  tubes  was  used. 

In  figure  2,  the  two  scintillator-photomultiplier  combinations 

22  60 

to  be  tested  are  separated  by  180°  centered  on  either  a  Na  or  Co 

22 

gamma  radiation  source,  as  the  radiation  from  the  Na  is  511  keV 

annihilation  radiation  emitted  at  180°  and  the  60Co  radiation  is 

1.173  and  1.332  MeV  cascade  radiation  emitted  in  a  distribution  peaked 

about  180°.  The  output  signal  from  one  of  the  detector  combinations 

23 

feeds  a  Constant  Fraction  of  Pulse  Height  Trigger  (CFPHT)  )  (to  give 
a  timing  signal  independent  of  pulse  height)  and  then  the  start  channel 
of  a  Time  to  Amplitude  Converter  (TAC)  while  the  other  signal  goes  into 


*Nash  &  Thompson  Ltd.,  Hook  Rise  South,  Tolworth,  Surbiton,  Surrey, 
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a  CFPHT,  is  delayed  a  set  amount  of  time,  and  then  goes  into  the  TAC 
stop  channel.  Pulse  height  requirements  are  provided  by  the  slow  out¬ 
puts  of  the  detectors.  These  signals  are  amplified  and  fed  into  Single 
Channel  Analysers  (SCA)  which  pass  only  signals  of  a  certain  size  range. 
Those  signals  which  are  not  blocked  by  the  SCA  are  then  used  to  strobe 
the  output  of  the  TAC.  The  ADC  presents  a  time  spectrum  which  indicates 
the  performance  of  the  detector  system  by  the  shape  and  the  width  of  the 
gamma  coincidence  peak  produced.  To  improve  the  resolution  and  extend 
the  applicable  time  region  to  shorter  lifetimes  an  effort  was  made  to 

obtain  better  photomultiplier  tubes  and  scintillators.  A  new  fast  scin- 

1  4 

tillation  material  1  -4  -bis  (2-butyloctyloxy)-P-quaterphenyl  as  reported 

24  * ** 

in  a  paper  by  Lynch  )  was  ordered  from  E.  Merck  and  was  dissolved  in 
deoxygenated  toluene  and  encapsulated  in  1  1/2"  diameter  cylinders,  1/2, 
1,  and  11/2  inches  high  by  Nuclear  Enterprises  Ltd  in  the  hope  that 
improved  timing  characteristics  might  be  achieved.  Two  new  types  of 
photomultiplier  tubes  were  also  tested, a  RuLlips  XP1210  and  an  improved 
version  of  the  RAC  8575  tube,  the  C31000D. 

The  first  test  conducted  consisted  of  mounting  1"  x  1"  diameter 
Naton  scintillators  on  two  new  C31000D  tubes.  It  was  expected  that  the 
peak  obtained  with  these  tubes  would  have  been  narrower  than  the  peaks 
obtainable  with  the  8575' s  due  to  the  much  improved  secondary  emission 

*E.  Merck  AG,  61  Darmstadt  2,  Postfach  4119,  Germany 

**Nuclear  Enterprises  Inc.,  935  Terminal  Way,  San  Carlos,  California, 
94070,  U.  S.  A. 
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60 

ratio  of  the  first  dynode  of  the  C31000D,  but  the  Co  FWHM  was  only 
210  psec  at  best,  about  the  same  figure  attainable  with  the  8575 ’s. 

A  Phillips  XP1210  photomultiplier  tube  had  been  obtained  from 
the  manufacturer  and  the  published  risetime  figures  are  better  than 

those  of  the  C31000D  (C31000D  10%  to  90%  anode  pulse  risetime  is  2.1  x 

-9  -9 

10  sec,  XP1210  risetime  is  1.2  x  10  sec).  This  tube  was  expected 

to  improve  the  resolution.  However,  when  one  of  the  C31000D's  was 

replaced  by  the  XP1210  the  FWHM  increased  to  ~270  psec.  Since  the 

overall  gains  of  the  C31000D  and  XP1210  are  about  the  same  (10^) ,  a 

possible  explanation  for  the  absence  of  any  timing  improvement  with  the 

latter  is  that  its  photocathode  efficiency  is  lower  than  the  efficiency 

of  the  C31000D  photocathode,  i.e.  the  luminous  sensitivity  of  the 

C31000D  is  85  yA/lm  measured  with  a  tungsten  filament  lamp  at  a  color 

temperature  of  2870°  K  and  the  XP1210  has  a  luminous  sensitivity  of 

45  yA/lm  from  a  similar  lamp  at  a  color  temperature  of  2854°  K.  The 

spectral  response  of  the  XP1210  is  Sll  while  the  C31000D  spectral 

response  is  very  close  to  S8*.  The  corrected  ratio  of  the  luminous 

sensitivities  and  hence  photocathode  efficiencies  of  the  C31000D  to  those 

of  the  XP1210  is  1.9. 

In  order  to  understand  the  probable  cause  of  the  lack  of  timing 
improvement  with  the  C31000D  over  the  8575,  the  operation  of  the  photo¬ 
multiplier  must  be  briefly  noted.  Consider  a  short  pulse  of  light 
striking  the  photocathode.  Electrons  are  emitted  in  a  time  distribution 
peaked  about  a  central  time  value  and  these  are  accelerated  and 
strike  the  first  dynode  causing  multiple  secondary  emission.  These 
electrons  then  strike  the  second  dynode,  again  causing  multiple  secondary 


*RCA  Tube  Handbook  HB-3  Vol  3,  1950. 


p- 

.no/  Jul OB91  *  3  •avoiCAoJ:  o3 

*  3'  •  Si  t  o  h b*j  £■  o  £>  cf  '  T-r;  =>■/:?  ,  J3  ild  o:-  i»bio  "i 


28 


emission.  This  process  is  repeated  along  the  dynode  chain  until  a  large 
output  pulse  is  produced  at  the  last  dynode  and  the  anode  of  the  tube. 

For  optimum  system  resolution  the  time  that  should  remain  as 
constant  as  possible  is  the  time  between  the  gamma  ray  striking  the 
scintillator  and  the  photomultiplier  tube  output  pulse  peak  time.  The 
timing  properties  of  the  tube  are  highly  dependent  upon  the  number  of 
photoelectrons  produced  in  the  first  one  or  two  stages  of  the  tube  in 
the  following  way.  If  many  electrons  are  produced  from  the  photocathode 
they  will  form  a  smooth  curve  in  a  time  distribution  with  a  well  defined 
peak  position.  While  the  number  of  electrons  arriving  in  a  portion  of 
this  curve  of  time  At  may  vary  slightly  from  pulse  to  pulse,  because  of 
the  large  number  of  electrons,  the  curve  shape  stays  pretty  much  the 
same,  and  the  peak  position  (central  time  value)  remains  the  same  from 
pulse  to  pulse. 

During  the  multiplication  stages,  each  photoelectron  will  have 
its  "descendants"  spread  in  their  own  time  distribution  curves,  but 
summed  at  the  output  pulse  they  still  form  the  smooth  curve  with  a 
constant  peak  timing  from  pulse  to  pulse, 

However,  with  few  electrons  emitted  from  the  photocathode,  a 
smooth  curve  is  not  formed  and  the  peak  position  will  change  from  pulse 
to  pulse,  causing  a  substantial  timing  jitter  and  loss  of  resolution. 

From  the  above,  it  can  be  seen  that  if  the  first  stage  or  two 
of  the  tube  can  be  made  to  have  very  high  statistics  in  the  electrons, 
better  timing  results,  and  while  both  the  8575  and  C31000D’s  have  good 
efficiency  in  the  photocathodes,  the  latter  tube  endeavours  to  become  a 
more  precise  timing  instrument  by  using  a  high  gain  gallium-phosphide 
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first  dynode.  This  should  have  the  effect  of  presenting  the  second 
dynode  with  smoother  pulse  than  in  the  8575.  However,  since  the  tubes 
behaved  similarly,  it  is  assumed  that  this  effect  is  not  as  significant 
between  the  first  and  second  dynode  as  between  the  photocathode  and 
first  dynode.  This  is  thought  to  be  because  the  flight  path  (and  hence 
straggling)  is  much  greater  between  the  photocathode  and  the  first 
dynode  than  between  the  first  and  second  dynode. 

While  this  testing  was  going  on,  a  few  runs  were  made  to  deter¬ 
mine  the  optimum  operating  voltages  for  the  photomultiplier  tubes.  These 
will,  of  course,  vary  with  gamma  ray  energy,  scintillator  type  and  from 

tube  to  tube,  but  for  the  1"  x  1"  diameter  Naton  scintillator  and  a 
60 

Co  source  the  tube  available  seemed  to  perform  optimally  at  about 
2300  V  for  the  8575’ s  and  about  2050  V  for  the  C31000D’s. 

At  this  point  in  the  bench  testing  procedure  the  quaterphenyl 
scintillators  arrived  from  the  manufacturer  and  testing  of  them  began. 
The  two  smallest  detectors  were  first  mounted,  the  1/2"  x  1  1/2"  dia¬ 
meter  on  an  8575  and  the  1"  x  1  1/2"  on  a  C31000D  photomultiplier  tube. 

The  quaterphenyl  scintillators  were  expected  to  reduce  the  FWHM  apprec- 

24 

iably  as  Lynch’s  paper  )  had  given  light  decay  time  constants  for  them, 

which  were  substantially  less  than  those  for  Naton  (60%  less  for  the 

energy  transfer  mean  time  and  20%  less  for  the  light  emitting  decay 

mean  lifetime).  However,  the  best  FWHM  which  could  be  obtained  was 

60 

207  psec  from  Co  (figure  3).  To  find  the  reason  for  the  lack  of 
improvement  in  the  FWHM,  a  pulse  height  study  was  carried  out  by  mount¬ 
ing  a  Nal  crystal,  a  Naton  crystal,  and  a  quaterphenyl  scintillator 
successively  on  the  same  photomultiplier  tube  (8575)  and  observing  the 
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pulse  height  by  recording  a  Na  Compton  spectrum  from  the  slow  channel. 
A  plot  of  these  pulse  height  spectra  is  given  in  figure  4.  Lynch’s 
paper  states  that  the  pulse  strengths  should  be  in  the  ratio  of  9. 0/1. 4/ 
2.1  for  NaI/Naton/l\4^-bis  (2-butyloctyloxy)-P-quaterphenyl .  The  data 
of  figure  4  gives  9.0/1.78/1.39  for  the  same  ratios.  The  scintillators 
were  constructed  according  to  Lynch’s  formula  and  as  the  above  results 
should  have  shown  characteristics  similar  to  his  published  figures,  it 
is  surmised  that  the  deoxygenation  of  the  scintillators  was  probably  not 
sufficient  during  the  manufacturing  stage.  The  resultant  reduced  pulse 
height  from  the  quaterphenyl  may  be  the  reason  for  the  FWHM's  from  these 
scintillators  not  being  smaller  than  those  from  Naton  crystals. 

As  a  further  investigation,  the  pulse  heights  of  the  three 
quaterphenyl  scintillators  ((1/2",  1",  2")  x  1  1/2"  diameter)  were  com¬ 
pared  by  successively  mounting  the  scintillators  on  the  same  8575  as  in 
the  Nal  v.s.  Naton  v.s  quaterphenyl  test.  The  lower  Compton  edge 
results  are  plotted  in  figure  5.  From  the  plots  it  was  noted  how  the 
size  of  the  scintillator  affected  the  pulse  height  and  this  suggested 
a  further  test.  Six  small  1"  diameter  cylinders  were  cut  from  the  same 
block  of  Naton,  two  1"  high,  two  1/2"  high  and  two  1/4"  high.  These 
were  highly  polished  and  mounted  in  MgO  lined  plastic  holders.  (It 
should  be  noted  that  polished  Naton  surfaces  should  not  be  touched  with 
any  cleaning  substance  except  absorbent  cotton  and  plastic  cleaner  and 
even  then  very  lightly  and  sparingly.)  The  scintillators  were  mounted 
in  pairs  of  the  same  size  on  RCA  8575  photomultiplier  tubes  and  the 
circuit  of  figure  2  was  again  used.  The  1"  x  1"  diameter  Naton  scin¬ 
tillators  were  first  tried  and  with  2200  V  on  each  tube  and  20%  windows 
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on  the  SCA  s  an  optimum  FWHM  of  193.6  +8.1  psec  was  obtained.  This 
is  plotted  in  figure  6. a.  (These  FWHM  were  fitted  using  a  new  version 
of  the  program  Tepel  provided  by  Mr.  D.  Orachesky.)  With  the  same 
window  size  the  1/2"  x  1"  diameter  Naton  scintillators  were  mounted, 
and  following  a  careful  calibration  carried  out  by  inserting  known 
delays  into  the  stop  signal  line  and  using  peaks  thus  created  to 
calibrate  the  ADC,  FWHM  of  165.0  +6.9  psec  was  achieved.  This  is  plotted 
in  figure  6.b.  The  1/2"  thick  Naton  cylinders  were  removed  and  replaced 
by  1/4"  thick  pieces  and  with  the  windows  reset  to  20%,  a  FWHM  of  162.9 
+6.8  psec  was  achieved.  This  is  plotted  in  figure  6.c.  The  measure¬ 
ment  with  the  1/2"  x  1"  diameter  Naton  scintillator  was  repeated  a  week 
later  as  a  check  and  at  that  time  FWHM  under  the  same  conditions  was 
170.4  +5.1  psec.  Unfortunately,  time  did  not  permit  further  study  of 
the  variation  of  peak  width  with  scintillator  size  at  this  time. 

Bench  testing  was  stopped  to  prepare  for  a  lifetime  measurement 
45 

run  on  a  Sc  target  using  the  Van  de  Graaff  accelerator.  The  detection 

system  chosen  on  the  basis  of  the  data  of  the  bench  testing  was  a  1"  x 

1  1/2"  diameter  quaterphenyl  detector  mounted  on  a  8575  photomultiplier 

tube.  Initially,  it  was  placed  16  cm.  from  the  target,  but  was  moved 

back  to  95  cm.  from  the  target  to  reduce  the  high  counting  rate,  when 
45 

a  very  thick  Sc  target  was  inserted  midway  through  the  run.  This 
scintillator  was  chosen  as  it  best  combined  fast  timing  characteristics 
with  a  large  active  volume  to  maximize  the  efficiency  of  the  experiment. 
This  detection  system  was  used  in  the  TWOD  data  collection  mode  des¬ 
cribed  above.  For  this  experiment  (and  others)  Dr.  W.K.  Dawson  of  the 
University  of  Alberta  has  provided  a  flexible  TWOD  program  with  variable 
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bin  size  and  number  which  ran  on  the  laboratory  Honeywell  DDP-516 

computer  using  the  display  unit  of  the  laboratory  SDS  920  computer. 

In  this  experiment,  sixteen  bins  of  256  channels  each  were  used  with 

the  time  calibration  of  74.85  psec/channel . 

This  time  calibration  was  done  by  a  method  basically  developed 
25 

by  Hatcher  )  and  modified  by  the  laboratory  staff  at  the  University 
of  Alberta.  The  output  of  a  10  MHz  sine  wave  generator  is  fed  into  a 
fast  zero  crossing  discriminator  to  give  a  train  of  pulses  100  nsec 
apart.  Two  outputs  from  the  crossing  discriminator  are  used,  one  to 
feed  the  start  of  the  Time  to  Amplitude  Converter  and  the  other  to  feed 
the  stop  after  passing  through  a  gating  device.  The  gate  is  controlled 

by  the  output  of  a  pulse  train  generator  whose  rate  is  adjusted  so  that 
the  output  of  the  TAC  will  be  pulses  of  amplitude  corresponding  to  a 
random  multiple  of  100  nsec.  When  collected  by  the  computer  and  dis¬ 
played,  the  time  per  channel  of  the  ADCA  spectrum  could  be  ascertained. 

A  problem  arose  in  producing  finer  calibrations  as  at  least  two  peaks 
must  show  on  the  display  for  a  calibration  to  be  made.  With  the  time 
scale  that  the  TAC  was  to  be  run  on  this  could  not  be  achieved,  so 
several  delays  were  calibrated  using  a  suitable  TAC  scale  and  these  were 
used  to  calibrate  the  correct  TAC  scale  noting  the  shift  in  channels  of 
a  single  artificially  generated  peak  with  the  insertion  of  the  delay 
into  the  stop  channel. 

The  calibration  of  the  bin  energy  (essentially  the  calibration 

of  the  ADCB  spectrum)  was  done  by  recording  the  Compton  spectrum  from 
22 

a  Na  source.  Assuming  linearity,  each  bin  covered  approximately  250 


keV  of  the  energy  spectrum. 
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For  each  run  an  ungated  Ge(Li)  spectrum  of  the  gamma  radiation 
was  also  recorded  with  a  30  cc.  active  volume  Ge(Li)  detector*  placed 
30.4  cm.  from  the  target.  After  taking  data  collection  runs  at  increas¬ 
ing  bombarding  energy  steps,  a  comparison  of  the  Ge(Li)  gamma  spectrum 
taken  at  a  given  energy  with  the  spectrum  taken  at  the  immediately  pre¬ 
ceding  energy  step  helped  to  identify  gamma  rays  by  defining  the  energy 
at  which  they  are  first  produced  and  hence  their  nuclear  level  of  origin. 
This  technique  is  of  particular  use  when  dealing  with  a  reaction  with 
a  negative  Q  value,  i.e.  a  reaction  which  requires  a  positive  energy 

input  if  it  is  to  go  at  all.  The  reaction  studied  in  this  experiment, 

45  45  * 

Sc(p,n)  Ti  is  such  a  reaction  with  a  Q  value  of  -2.841  MeV.  This 

reaction,  suggested  as  a  candidate  for  such  a  study  by  Dr.  K.V.K. 

45 

Iyengar,  involves  bombarding  Sc  with  protons  and  ignoring  the  emitted 

45 

neutrons  while  observing  the  decay  of  the  excited  states  of  Ti  which 

45  45  * 

are  formed.  The  competing  reaction  Sc(p,p)  Sc  cannot  be  prevented 

and  must  be  accounted  for  in  all  results  such  as  lifetime  and  Ge(Li) 

45  42 

spectra.  Fortunately,  the  competing  reaction  Sc(p,a)  Ca  was  not  too 

42 

prominent  a  reaction  due  to  the  high  Coulomb  barrier  which  the  Ca 
presents  to  the  a  cluster.  Other  competing  reactions  from  protons  on 
^Sc  require  input  energies  in  excess  of  10  MeV  and  so  are  not  consider¬ 
ed  in  this  work. 

The  targets  were  generally  the  so  called  "goop"  or  "gloop" 
targets  which  consist  of  a  mixture  of  polyurethane  dissolved  in  benzene 
into  which  the  finely  powdered  target  material  is  mixed.  The  matrix 
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holding  the  target  material  is  thus  composed  of  only  carbon,  oxygen  and 
hydrogen,  none  of  which  have  levels  with  lifetimes  that  would  interfere 
with  any  experiments  being  carried  out.  This  brings  up  an  important 
experimental  point,  that  the  beam  should  not  be  able  to  hit  anything 
except  the  target  or  pure  gold  as  other  materials  might  present  "long 
lived"  gamma  rays  (e.g.  tantalum  10.6  ns,  482  keV  level)  which  would 
interfere  with  the  gammas  of  interest.  The  target  materials  were  thus 
"glooped"  onto  0.020"  gold  backings  which  were  in  turn  mounted  on  a 
2  1/2"  diameter  disc  (with  sping  clips  on  one  face  to  hold  the  gold- 
backed  target) .  This  disc  fitted  over  the  end  of  the  beam  tube  and 
formed  a  butt  vacuum  seal  with  that  end  and  was  also  secured  in  place 
by  a  plastic  retaining  ring.  Unfortunately,  the  disc  was  constructed 
of  tantalum  as  it  was  designed  for  another  experiment  which  did  not 
involve  lifetime  studies  and  hence  did  not  require  this  same  precaution. 
To  solve  this  problem  a  gold  annular  shield  was  constructed  to  fit 
inside  the  beam  tube  completely  covering  the  interior  end  except  for  a 
1/2"  diameter  opening  exposing  the  target  to  the  beam. 

A  stub  extending  from  the  tantalum  disc  through  a  slot  on  the 
retaining  ring’s  inner  face  formed  a  connection  for  the  Faraday  cup  and 
cooling  of  the  targets  was  provided  by  air  blown  against  the  back  of  the 
disc.  On  occasion  when  a  thick  metal  target  (which  would  not  let  the 
beam  strike  the  tantalum)  was  used,  the  target  was  mounted  without  gold 
backing  directly  on  the  disc  to  provide  better  cooling. 

The  run  was  begun  with  a  10B  target  instead  of  ^5Sc  to  check 
the  operation  of  the  circuit.  The  data  of  figure  7  (a  natural  logarithm 
plot)  was  taken  with  a  charge  collection  of  3972  uc  in  just  over  forty- 
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five  minutes.  The  lifetime  tails  in  bins  0,  1  and  2  are  due  to  gamma 

radiation  emitted  by  the  717  keV  level  of  10B (p ,p ' ) 10B*  reaction.  These 

curves  were  analysed  on  the  University  of  Alberta  IBM  360  computer  by 

using  a  weighted  least  squares  calculation  due  to  Ferguson^) .  The  work 

was  done  on  the  Nuclear  Research  Centre's  remote  computer  terminal  using 

the  APL  language  and  copies  of  the  fitting  and  display  programs  are 

included  in  Appendix  III.  The  analysis  gives  a  mean  life  of  1026.3  + 

8.6  psec  for  the  mean  nuclear  lifetime,  T  .  of  the  717  KeV  level  of  10B. 

m 

This  corresponds  to  a  mean  half  life  of  711.6  +6.0  psec  which  agrees 

1  Q 

favorably  with  the  published  figure  of  0.71  nsec  ). 

After  the  run  was  completed  a  thin  gold-backed  ^Sc  target 

was  inserted.  (It  should  be  noted  that  changing  targets  with  this 

apparatus  could  be  accomplished  very  quickly  as  time  was  not  required  to 

depressurize  a  large  vacuum  vessel.)  In  order  to  carry  out  a  lifetime 

45 

survey  of  the  Ti  levels  it  was  necessary  to  have  a  level  scheme  and 

45 

this  was  provided  by  Dr.  K.V.K.  Iyengar  from  his  work  on  the  Sc(p,n) 

26 

reactions.  These  levels  are  shown  on  figure  8  ).  With  the  -2.841  MeV 

Q  value  proton  bombarding  energies  of  3.15,  3.35,  3.955,  4.25,  4.3,  4.45 
and  4.55  MeV  were  chosen,  each  energy  being  just  above  the  threshold  for 
the  0.40,  0.330,  0.740,  1.225,  1.350,  1.465  and  1.515  MeV  levels  respec¬ 
tively.  The  data  is  not  presented  in  order  of  ascending  proton  bombard¬ 
ing  energy  but  rather  in  chronological  order  as  this  makes  evident  the 
unsteadiness  and  general  decrease  of  the  beam  resolution  with  time.  The 
beam  was  initially  adjusted  for  optimum  pulse  size  and  sharpness  and  at 
the  optimum  resolution,  the  shortest  lifetime  which  could  be  measured  was 
92.5  +  5.1  psec  (this  being  the  value  obtained  by  applying  equation  2.b 
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to  the  left  hand  side  of  the  prompt  peak) .  The  peaks  in  the  upper  bins 
of  figure  9a  (E  =  3.955  MeV)  were  taken  as  the  prompt  peaks  and  appli¬ 
cation  of  the  least  square  fitting  program  gave  their  inverse  slopes 
(lifetime)  as  in  table  III. 


TABLE  III 


Bin 

11 

94.2 

+ 

5.2 

psec 

Bin 

12 

107.7 

+ 

5.6 

psec 

Bin 

13 

92.5 

+ 

5.1 

psec 

Bin 

14 

98.2 

+ 

6.2 

psec 

After  three  runs  had  been  taken  with  the  thin  target  (E^  =  3.955,  4.3, 
4.45  MeV;  figures  9a,  9b,  and  9c  respectively),  a  very  thick  (0.005") 
scandium  foil  was  inserted.  The  detector  was  moved  back  and  another  run 
was  taken  at  4.45  MeV,  shown  in  figure  9d.  The  data  taken  at  E^  =  4.55 

MeV  is  shown  in  figure  9e,  4.25  MeV  in  figure  9f,  3.35  MeV  in  figure  9g 

and  3.15  MeV  in  figure  9h.  The  actual  data  taking  was  comparatively 

uneventful  apart  from  slight  trouble  with  the  beam  conditions  drifting 
(cf.  figure  9b). 

Due  to  this  beam  drifting  the  prompt  peaks  taken  initially  could 
not  be  used  throughout  the  run  and  so  analysis  of  the  runs  could  not 
proceed  by  the  centroid  shift  or  by  the  conventional  method  given  in 
section  III.  As  the  prompt  peak  could  not  be  subtracted,  the  shortest 
lifetimes  are  uncertain  as  they  could  not  be  separated  from  the  prompt 
peak  contribution.  The  method  of  analysis  used  consisted  of  visually 


scanning  the  natural  logarithmic  plots  of  the  data  to  find  any  peaks 
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whose  left  hand  slopes  were  less  steep  than  their  right  hand  slopes 
indicating  the  presence  of  a  measureable  lifetime.  Such  slopes  existed 
with  all  the  runs  in  at  least  bins  0  to  6 .  Since  most  of  these  bins 
also  had  noticeable  longer  lifetimes  as  well,  the  unfolding  procedures 
and  programs  described  above  were  employed.  These  gave  the  results 
shown  in  table  IV. 


TABLE  IV*  ** 


Run;  E 

P 


Bins 


T  (psec) 
m 


4/6/70/3;  4.45  MeV 


1 

7500 

± 

1200 

470 

+ 

100 

104.7 

+ 

0.9 

1 4  - 

15 

2800 

+ 

780 

420 

+ 

160 

97.4 

+ 

1.6 

4/6/70/4; 

4.45  MeV 

- 

1 

6100 

+ 

1300 

594 

+ 

27 

166.1 

+ 

0.6 

2 

1265 

+ 

82 

406 

+ 

41 

154.8 

+ 

0.8 

3 

690 

+ 

46 

170.6 

+ 

1.7 

4 

920 

+ 

150 

174.2 

+ 

1.7 

5 

707 

+ 

56 

158.8 

+ 

2.1 

6 

750 

+ 

110 

182.7 

+ 

6.0 

7 

760 

+ 

190 

147.4 

+ 

3.4 

8 

621 

+ 

73 

152.9 

+ 

3.8 

11 

600 

+ 

130 

140.6 

+ 

5.1 

14 

231 

+ 

17 

117 

+ 

11 

*Bin  numbers  are  assigned  with  the  lowest  bin  designated  bin  0. 

JL  JL 

Errors  include  only  statistical  errors. 


4/6/70/5;  4.55  MeV 


2 

5100 

+ 

980 

482 

+ 

42 

180.9 

+ 

1.0 

3 

8800 

+ 

3700 

556 

+ 

59 

198.2 

+ 

1.3 

5 

656 

+ 

52 

191.1 

+ 

1.5 

6 

776 

Hr 

99 

195.8 

+ 

3.5 

7 

647 

+ 

84 

181.8 

+ 

4.1 

8 

607 

+ 

85 

186.0 

+ 

4.3 

11 

475 

+ 

69 

186.6 

+ 

6.0 

13 

469 

+ 

85 

161.4 

+ 

7.7 

14 

279 

+ 

22 

164.7 

+ 

2.0 

4/6/70/6; 

4.25 

MeV 

2 

3800 

+ 

1000 

442 

+ 

33 

153.8 

+ 

1.2 

3 

3800 

+ 

1400 

468 

+ 

47 

161.8 

+ 

1 . 6 

4 

4400 

+ 

1700 

467 

+ 

75 

161.7 

+ 

3.1 

5 

520 

+ 

29 

158.0 

+ 

3.5 

6 

744 

+ 

98 

166.1 

+ 

4.8 

11 

930 

+ 

250 

163.6 

+ 

5 . 6 

4/6/70/7; 

3.35 

MeV 

2 

4000 

+ 

1300 

573 

+ 

75 

134.5 

+ 

1.7 

3 

2460 

+ 

980 

455 

+ 

92 

110.7 

+ 
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The  longest  lifetimes,  those  in  the  2500  to  6000  psec  range,  are 
due  to  attempts  to  fit  straight  lines  to  the  very  longest  lifetime 
tails  and  as  the  statistics  are  very  poor  in  that  region,  the  errors 
are  correspondingly  large.  Thus,  with  the  longest  lifetimes  uncertain 
and  the  shortest  ones  being  interfered  with  by  the  still-present  prompt 
peak,  only  one  lifetime  can  be  inferred  from  the  data.  This  is  the 
565  +  20  psec  lifetime  found  in  all  the  runs  and  bins  up  to  bin  13 
inclusive.  This  gamma  ray  is  between  2.8  and  3.3  MeV  in  energy  accord- 

fpref,  I'ntTfly 

ing  to  the  ADCB  calibration  and  the  position  of  its  Compton  edge  (bin 
II 

It). 

The  origin  of  this  gamma  ray  is  now  considered.  It  was  estab¬ 
lished  above  that  the  only  contributing  reactions  are  ^^Sc(p,p'  )^Sc* 

45  45  * 

and  Sc(p,n)  Ti  .  If  the  reaction  exciting  the  565  psec  lifetime 
level  were  the  (p,n)  reaction,  the  gamma  ray  in  question  would  not  be 
seen  at  all  if  the  proton  bombarding  energy  were  below  the  threshold  for 
the  creation  of  that  particular  state  (cf.  figure  8).  Above  threshold, 
the  gamma  ray  would  be  seen  at  every  bombarding  energy  due  to  the  , 

continuum  of  proton  energy  caused  by  straggling  in  the  target  material; 
and  its  appearance  as  the  bombarding  energy  crossed  the  threshold  would 
be  quite  dramatic.  As  the  "565  psec  gamma  ray"  is  evident  in  run  4/6/ 
70/8  (Ep  »  3.15  MeV)  which  is  just  above  threshold  for  the  40  keV  level, 
this  would  imply  that  the  40  keV  level  has  a  565  psec  mean  lifetime  if, 
in  fact,  the  (p,n)  reaction  were  responsible  for  exciting  the  lifetime 
state.  However,  the  gamma  ray  in  question  is  approximately  3  MeV  in 
energy  as  mentioned  above.  This  implies  that  the  Sc(p,n)  Ti*  reac¬ 
tion  did  not  feed  the  level  having  this  measureable  mean  lifetime  and 
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it  is  therefore  concluded  that  one  of  the  levels  of  Sc  between  2.8 
and  3.2  MeV  has  a  lifetime  of  565  +  20  psec.  The  Ge(Li)  detectors  used 
for  these  runs  had  very  low  efficiencies  above  2.5  MeV  and  their  spectra 
could  not  be  used  as  an  aid  to  further  identification.  As  a  result,  to 
determine  precisely  which  level  has  the  measureable  lifetime  will 
require  a  further  study  perhaps  aided  by  the  refinements  to  be  discussed 
in  section  V. 

From  the  above  discussion,  it  may  also  be  concluded  that  the 
45 

Ti  has  no  strongly  (p,n)  excited  levels  with  mean  lifetimes  longer 
than  180  psec  (from  run/bin  4/6/70/5/2)  up  to  the  1.52  MeV  level  of 


excitation. 
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V.  CONCLUSIONS 


The  data  and  arguments  of  section  IV  give  an  upper  limit  of 
180  psec  on  the  mean  lifetimes  of  any  strongly  excited  levels  of 
^5Ti  up  to  1.52  MeV  excitation  energy  excited  by  the  ^5Sc(p,n)^5Ti* 
reaction.  In  addition,  a  new  lifetime  of  565  +  20  psec  has  been 
assigned  to  the  excitation  region  of  ^5Sc  between  2.8  and  3.2  MeV 
and  modifications  of  the  above  technique  to  positively  identify  the 
responsible  level  are  discussed  below. 

27 

While  no  evidence  of  the  previously  published  )  1100  psec 

45 

lifetime  level  at  376  keV  in  Sc  was  found,  it  was  carefully  looked 
for.  However,  the  analysis  showed  that  the  tail  of  the  565  psec  gamma 
ray  merging  into  the  background  fully  masked  the  region  where  this 
lifetime  tail  should  have  been  found  in  the  lower  bins.  Probably  the 
intensity  of  the  1100  psec  gamma  ray  is  quite  weak.  Better  statistics 
may  show  a  separation  of  the  two  lifetimes  if  the  longer  one  is  a  real 
feature.  Along  this  line,  general  improvements  to  the  technique  and 
the  analysis  are  to  be  discussed. 

The  major  problem  with  the  analysis  was  the  separation  and 
identification  of  the  various  components  in  the  lifetime  tails.  This 
was  due  to  the  poor  statistics  of  the  data  and  the  lack  of  a  fine 
enough  energy  calibration  in  the  upper  bins.  The  obvious  solution  to 
the  first  problem  is  to  extend  data  collection  time  and  to  use  thick 
targets  throughout  the  experiment  (care  must  be  taken  with  targets 
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that  become  substantially  strong  long-lived  radiation  sources).  However, 
as  the  collection  times  were  in  the  order  of  thirty  minutes  for  a  run 
and  the  number  of  counts  in  the  lifetime  tails  should  be  increased 
manifold,  the  running  time  required  might  become  prohibitive.  While 
long  runs  are  out  of  necessity  being  considered,  it  might  be  well  worth¬ 
while  to  give  consideration  to  particle-gamma  coincidence.  This  would 
positively  identify  a  gamma  ray  by  gating  the  collected  timing  peaks 
with  coincident  emitted  particles  of  definite  energy.  With  a  knowledge 
of  the  level  scheme  from  gamma  ray  spectroscopy,  choosing  only  those 
gammas  from  a  level  defined  by  a  given  outgoing  particle  energy  would 
allow  a  lifetime  study  to  be  made  level  by  level.  In  addition,  without 
competing  gamma  rays  from  other  levels  and  reactions,  the  numbers  of 
counts  in  the  lifetime  tails  required  in  order  to  obtain  a  value  for 
a  mean  lifetime  would  be  less  than  the  numbers  collected  for  this  exper¬ 
iment,  small  as  they  were.  , 

The  second  problem  may  be  overcome  with  an  accurate  energy 

12 

calibration  in  the  higher  bins  using  the  C  first  excited  state  at 
4.43  MeV  by  the  1^C(p,p ' )^2C*  reaction. 

Another  problem  was  noted  during  analysis.  While  the  particle 
gamma  coincidence  technique  mentioned  above  would  eliminate  the  need 
for  a  prompt  peak,  when  the  noncoincidence  technique  is  used  as  in  this 
work,  a  method  for  removal  of  the  prompt  peak  is  found  to  be  necessary 
to  investigate  the  shorter  lifetimes.  From  observations  while  running, 
drifts  in  the  system  which  affect  the  resolution  appear  to  occur  at  a 
comparatively  slow  rate  so  that  one  suggestion  for  obtaining  a  prompt 
peak  is  to  place  the  target  material  on  half  the  surface  of  a  rapidly 
rotating  cylinder,  the  other  half  of  which  is  to  be  covered  by  a 
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material  producing  only  prompt  gamma  radiation.  With  the  lifetimes 

being  considered  in  this  work,  the  data  from  the  detector  could  be 
gated  into  the  data  collection  region  during  the  exposure  of  the 
"prompt"  material  to  the  beam  and  into  a  second  region  for  the  actual 
data  collection  while  the  target  material  was  exposed.  This  would 
allow  a  nearly  simultaneous  collection  of  both  the  data  and  a  prompt 
peak  under  the  same  beam  conditions.  Other  methods,  based  on  deriving 
prompt  peaks  simultaneously  with  the  data,  are  currently  being  explored 
at  the  Nuclear  Research  Centre. 

Perhaps  the  greatest  aid  to  this  work  would  be  an  on-line  anal¬ 
ysis  program  utilizing  either  the  Honeywell  DDP-516  computer  or  the  SDS- 
920.  This  type  of  program  would  allow  an  evaluation  of  the  statistics 
to  be  taken  during  such  an  experiment  and  on  the  basis  of  the  evaluation, 
the  experiment  would  be  extended,  modified  or  even  terminated  (at  a 
tremendous  saving  in  time  and  effort  to  the  experimentor ) . 

If  on-line  analysis  cannot  be  provided  (a  situation  which  would 
result  in  the  ambiguities  present  in  this  work)  it  is  suggested  that  a 
complete  analysis  program  working  from  data  on  magnetic  tape  be  devel¬ 
oped  with  provisions  for  prompt  peak  subtraction  and  multiple  lifetime 
extraction.  This  program  should  be  based  on  the  unfolding  alogrithm 
mentioned  in  section  IV  and  the  resultant  APL  programs  in  Appendix  III. 

In  conclusion,  it  is  noted  that  while  the  method  of  finding  and 
evaluating  nuclear  lifetimes  by  direct  electronic  measurement  does 
suffer  from  statistical  problems  as  performed  in  this  work,  the  above 
suggestions  should  eliminate  this  problem  and  make  the  method  much 
easier  to  carry  out  (and  more  practical  with  respect  to  the  time 
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required  for  data  analysis).  This,  combined  with  the  simplicity  of 
the  method  and  its  applicability  to  a  wide  range  of  target  nuclei 
should  make  it  a  useful  experimental  tool. 
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VI .  APPENDIX 


APPENDIX  IA>1) 


The  vector  potential  A  which  satisfies  the  wave-equation  in  free 


space 


[i.i] 


can  be  used  to  derive  expressions  for  an  electromagnetic  field  including 
the  electromagnetic  field  of  y-rays  emitted  by  excited  nuclei.  By  elim¬ 
inating  the  time  dependent  part  of  the  scalar  potential  the  electromag¬ 
netic  waves  are  seen  to  be  transverse  from 


div  A  =  0  . 


[1.2] 


The  electric  and  magnetic  fields  are  given  by 


[I. 3. a] 


E 


c  3t 


[I.3.b] 


H  =  curl  A 


In  quantizing  the  radiation  field  the  vector  potential  is  expanded 


as  a  linear  combination  of  the  simultaneous  eigenfunctions  of  equations 
1.1  and  1.2.  By  confining  the  radiation  to  a  specified  volume  and  im¬ 
posing  boundary  conditions  any  vector  potential  A  which  satisfies  the 
boundary  conditions  can  be  written  as  a  linear  combination  of  orthogonal 
waves . 


[1.4] 


7  -  - 
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Here  X  represents  the  various  quantum  numbers  associated  with  each  wave 
& 

(q*  anc*  ^x  defined  below) .  The  are  simultaneous  solutions  of 

the  vector  Helmholtz  equation 


+  kk 


=  o 


[1.5] 


(kx  being  the  eigenvalue  wave  number  characterizing  A  )  and  of  the  diver- 
gence  condition 

div  A ^  =  0  [ ! . 6 ] 

The  A^  are  normalized  for  convenience  to 

A^  ,  (r_)  •  A^  (_r) dv  =  4ttc26x,x  [1.7] 

V 


A  single  divergenceless  vector  Helmholtz  equation  can  be  derived 
by  combining  1.5  and  1.6  to  get 


curl  curl  A, 

A 


kk 


o 


[1.8] 


Solution  of  this  equation  will  be  discussed  below. 


The  field  amplitudes  q  and  q  arbitrarily  satisfy  the  harmonic 

A  A 

oscillator  equation 


dt 


2 


+ 


2 

03  q^ 


0 


[1.9] 


* 

where  a)  =  ck^  is  the  frequency.  The  q^  and  q^  represent  annihilation 
and  creation  operators  respectively  apart  from  a  multiplicative  time 
constant  and  a  harmonic  time  factor.  These  are  normalized  to  obey  the 
commutation  rules 


Vx'  ■  qx’qx 


0 


[I. 10. a] 
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[I.lO.b] 


*  *  *  * 

qxV  •  V'Ia  =  0 


"»v  -  ■  <£->*,., 

A 


[I.lO.c] 


The  matrix  elements  of  qA  and  q*  taken  between  an  ititial  state  with  n 

quanta  characterized  by  quantum  numbers  A  and  a  final  state  containing 
n£  quanta  also  characterized  by  A  are  given  by 


<nf ,  A  lqA  lni,A>  "  ni/h/2a)A1/2<Snf  ,n.  exP(-^xt) 


[I. 11. a] 


<nf ,AlqAlni,A>  "  (ni+1)*/2oj  1/26  exp(ia)  t) 

nf,ni+l  A 


[I.ll.b] 


The  radiation  field  energy  eigenvalues  are  now  found  in  the  fol¬ 
lowing  way.  The  above  expressions  are  used  to  change  the  radiation 
field  Hamiltonian 


into 


H  =  (8tt) 


-1 


i  2  i  i  2 

E  +  H  dV 


V 


H  =  l  2“a1*1x 


Thus,  from  I. 10. a  and  I.lO.b  we  get 


[1.12] 


[1.13] 


En  =  I  [1.14] 

A 

if  the  order  of  the  field  amplitudes  q^  and  q^  are  chosen  so  that  there 
is  zero  energy  when  no  quanta  are  present. 


In  the  following  work,  the  vector  potential  is  expanded  in  a 
series  of  standing  spherical  waves  centred  at  0.  A  sphere  of  radius 
is  chosen  as  the  volume  of  integration  with  the  boundary  conditions  that 
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the  tangential  components  of  A  vanish  at  the  surface  of  the  sphere.  The 
vector  potential  for  each  spherical  wave  is  constant  throughout  the  vol¬ 
ume  of  the  sphere  and  decreases  as  r  ^  for  large  r  from  the  centre.  Each 
of  the  waves  has  a  definite  angular  momentum  L  and  corresponds  to  the 
outgoing  spherical  wave  emitted  by  a  suitably  oscillating  2^  pole  at  0. 

A  box  quantization  with  plane  waves  is  not  used  as  the  plane  waves  do 
not  have  definite  angular  momentum  and  the  classical  multipole  oscillator 
correspondence  is  also  missing.  Before  finding  a  solution  to  the  diver¬ 
genceless  vector  Helmholtz  equation  (equation  II. 7),  solutions  to  the  scalar 
Helmholtz  equation 

Au  +  k^u  =0  [1.15] 

will  be  discussed.  It  can  be  shown  that 


u  =  i  (k  r)Y  (0td))  [1.16] 

with  integer  L  >  0  and  |m|  <  L,  satisfies  1.13.  The  1^(9, <j>)  are  the 
normalized  spherical  harmonics  and  the  JL(^Ar)  are  the  spherical  Bessel 
functions 


VkAr)  =  (2kT) 


TT  N  1/2 


JL+l/2(kAr) 


[1.171 


When  u(R  )  =  0  is  imposed  as  a  boundary  condition,  ^(k^R^  0 


which  requires 

k^  =  (n  +  (L/2))Rq 

where  n  is  a  positive  integer  for  Rq  >>  L7t/2k. 


[1.18] 


Since  solutions  to  the  vector  Helmholtz  equation  will  be  ex¬ 
pressed  in  terms  of  the  uLM’s  several  characteristics  of  the  latter 

will  now  be  noted. 

1)  The  u  form  a  complete  orthogonal  set  of  non-singular  solu 
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tions  to  the  scalar  Helmholtz  equation  with  the  orthogonality  conditi 
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[1.19] 


2)  Defining 


L.  =  -  x  grad. 


[1.20] 


the  satisfy 


L  u  =  Mu_.. 
z  LM  LM 


[I. 21. a] 


L  uLM  =  L(I*1)ulm 


[I . 21.b ] 


This  implies  that  the  u  are  eigenfunctions  of  the  orbital  angular 


momentum  operators 

3) 


u  (-r)  =  (-1)  u  (r) 
LM  LM  — 


[1.22] 


The  u^M  thus  have  even  (+1)  or  odd  (-1)  parity  depending  upon  the 
value  of  L.  The  solutions  to  the  divergenceless  vector  Helmholtz  equa¬ 
tion  can  be  expressed  in  terms  of  the  uT  ’s  as  mentioned  above,  but  for 

LM 

each  two  linearly  independent  solutions  are  chosen  where  a  is 

either  E  or  M  denoting  electric  or  magnetic  multipole  fields. 


=  kX1CLCUrl  — ULM 


[I. 23. a] 


aM 

A  =  tCLu 
— LM  L—  LM 


is  an  arbitrary  constant  and 


[1 . 23 .b  ] 
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The  ^LM  have  Pr°Perties  similar  to  the  u  .  A 

LM 

1)  The  form  a  complete  orthogonal  set  of  non-singular  solu¬ 
tions  to  the  divergenceless  vector  Helmholtz  equation.  The  orthogonal- 
ity  condition  is 
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L(Lfl) 

2k2 
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Equation  1.7,  the  normalization  condition  is  satisfied  if 


8ttu) 


2  1/2 


L(L+1)R 


2)  Defining  J  =  L  +  S 


and 


S  A  =  iz  x  A 

z—  — z  — 


[1.26] 


[1.27] 


[1.28] 


where  ^  is  a  unit  vector  in  the  z  direction  and  S  is  the  z  component 

z 

of  intrinsic  spin,  the  satisfy 


Km 


[I. 29. a] 


-  ^LM  L^L+1^M 


[1 . 29 .b  ] 


S2A^m=  S(S+1)A^m=  2A^m  [I.29.C] 

This  implies  that  the  elctromagnetic  field  associated  with  A^  con¬ 
tains  a  total  angular  momentum  !L  with  component  M  in  the  z  direction. 

differ  from  the  u^M  in  that  the  former  have  the  divi¬ 
sion  into  electric  and  magnetic  multipoles.  The  electric  and  magnetic 
fields,  E  and  H,  are  related  to  the  time  independent  vector  potentials  by 


3)  TheA^M 


aJ  no!  rbnoo  x3-* 

f  _ 
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[I. 30. a] 

— LH  =  CUrl 

[1 . 30 .b  ] 

From  this  it  follows  that 


=  =  iCL  CUrl  — ULM 

[1.31. a] 

=  "SlM  =  _kACL-ULM 

[1 . 31.b  ] 

which  satisfy  the  equation 

\A 
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APPENDIX  IIA’2) 


The  development  of  the  quantum  mechanical  transition  probability 
from  perturbation  theory  starts  with  the  total  Hamiltonian 


H=Hq+H'  [II. 1] 

being  the  sum  of  an  unperturbed  Hamiltonian  Hq  with  normalized  eigen¬ 
functions  u  and  eigenvalues  E 
n  n 


Hu  =  E  u  , 
on  n  n 


[II. 2] 


and  a  small  time  dependent  perturbation  H' .  Stationary  solutions  of  the 
Schrodinger  equation  do  not  exist  with  the  total  Hamiltonian  and  the 
time-dependent  equation 


if 1 


3t 


[II. 3] 


must  be  used.  In  order  to  solve  this  equation,  ip  is  expressed  as  an 
expansion  in  the  eigenfunctions  of  the  unperturbed  time-dependent  equa¬ 
tion  using  time  dependent  expansion  coefficients. 

-i E  t/fi 

ip  =  Y  a  (t)u  e  n  [11*4] 

Y  L  n  n 


where  £  is  to  be  taken  as  a  summation  over  the  discrete  set  of  eigen 
functions  and  integration  over  the  continuous  set.  Putting  II. 4  into 
II. 3  gives 


-tE  t/fi  -iE  t/h  -iE  t/h 

^anune  +IanEnV  "  -  l*n«0  +  ITL.S} 


Replacing  #QUn  with  E^u^,  left  multiplying  by  u^>  and  integrating  over 
all  space  gives 
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tfia^e 


CEL  t/fi  -i E  t/h 

k  r  n 

=  )a  e 


u,  tf'u  dx  . 
k  n 


[II. 6] 


For  convenience 


H'u  dx  =  H' 


V  un 


kn 


[II. 7] 


which  is  called  the  matrix  element  of  the  perturbation.  Further  defin¬ 


ing 


00. 


Ek  -  En 


kn  h 


[II. 8] 


equation  II. 6  becomes 


*-l  r  ^Wkn^ 

a.  =  (tfi)  1  I H!  a  e  kn 
k  ^  kn  n 


[II. 9] 


The  perturbation  approximation  consists  in  replacing  H'  by  \H' 
in  II. 1  and  II. 9  and  by  expressing  the  afs  as  a  power  series  in  X 


a 

n 


(0) 


n 


+  Xa 


(1) 


n 


X2a(2) 

n 


[II. 10] 


Equating  coefficients  of  equal  powers  of  X  and  setting  X  =  1  the  set  of 
equations  results 


a^0)  =  0;  a£s+1) 


l  H' 


f  \  t 

( s )  kn  _  rv  i  *) 

3.  0  }  S  - 

kn  n 


in. ii] 


The  first  equation  for  a^/s  shows  that  the  a^  ^  are  constant  in  time. 
Their  values  give  the  initial  conditions  of  the  problem  before  the  per¬ 
turbation  is  applied.  It  is  now  assumed  that  all  but  one  of  the  a^ 
are  zero,  implying  that  the  system  is  in  a  definite  unperturbed  energy 
state  at  the  time  the  perturbation  is  first  applied.  Thus 
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ak°>  =  6km  01  ak°>  =  6(k  '  m)  [IX. 12] 

depending  upon  the  discreteness  of  the  state  m.  Using  11.12  the  first 
order  equation  gives 


a^Ct)  =  (t*)"1 


•t  £(i0,  t ' 

^m(f)e  km  df 

— oo 


[11.13] 


When  #km  is  °f  finite  duration,  the  amplitude  of  a  state  u^  (k  =/  m) 
(after  the  perturbation  has  disappeared)  is  proportional  to  the  time 
Fourier  component  of  that  corresponds  to  u)^.  If  is  independent 

of  time  except  for  being  turned  on  at  time  0  and  off  at  time  t,  the 
first  order  amplitudes  at  the  time  t  is 


H,'  t'“knt  . 

(t) - Is- e  „  [11.14] 

Tl  CO. 

km 

Thus,  if  the  system  is  originally  in  state  m,  the  probability  of  finding 
the  system  in  state  k  at  time  t  is 


(1) 


From  the  plot  of  sin  ( 


(t)l2  - 

'  '"km1 

*2  2 
n  co. 

km 

0)  .  t 
km 

1  ?  2 

2o)  . 

km 

V .  S  .  CO. 

km 

km 


[11.15] 


A. 2 


(cf.  ))  it  can  be  seen  that 


k  (E^-Ejq)  is  proportional  to  t  if  i-s  approximately  independent  of  k. 


The  transition  probability  per  unit  time,  W,  is  found  by  assuming 

the  system  to  be  in  a  cubical  box  of  dimension  x  with  periodic  boundary 

conditions  at  the  walls.  The  eigenfunctions  un  now  form  a  discrete  set 

and  are  normalized  to  unity  in  the  volume  x3.  A  group  of  final  states  k 

with  E  ~  E  are  considered  and  it  is  assumed  that  varies  slowly  with 

k  m 


k.  The  density  of  final  states 


All 


dN 
dE  ” 


P  GO 


[11.16] 


is  assumed  to  a  slowly  varying  function  of  k  also.  Using  this,  the 
transition  probability  per  unit  time  to  some  state  k  of  the  final  group 
is 


W  = 


-1 


a^1)(t)|2p(k)dEk 


[11.17] 


where  the  integration  is  introduced  by  taking  x  large  enough.  The 
and  p  (k)  are  taken  outside  the  integral  as  they  are  slowly  varying  and 
most  of  the  contribution  comes  from  a  small  energy  range  about  E^,  giving 


4 \h: 

1  km 
fa 


roo 


p  00 

—  00 


sin^(o).  t/2) 
km 


GO. 


km 


[11.18] 


Since  go.  ~0,  the  integral  can  be  rewritten  as 
km 


f  00 


2 


-2.2,  TTt 
x  sin  xdx  =  — 


so  that 


w  =  f  p(k)KJ2 


or 


[11.19] 


[II. 20. a] 


[II . 20 .b  ] 


using  the  notation  of  equation  6,  section  II. 
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APPENDIX  IIIA*3) 
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VII.  GLOSSARY 


Analogue  to  Digital  Converter:  ADC,  an  electronic  device  which  assigns 
digital  information  to  an  input  analogue  signal  proportional  to 
the  input  signal  size. 

Compton  edge:  the  drop  off  in  counts  with  increasing  energy  in  an  energy 
spectrum  due  to  the  maximum  allowable  value  of  the  energy 
transfer  in  Compton  scattering 


T 

max 


hv 


1  + 


hV 


me 


A  point  two-thirds  of  the  way  down  the  Compton  edge  is  considered 
T 

max 

continuum:  the  energy  region  where  the  discrete  energy  states  of  a 

nucleus  are  so  close  together  that  their  widths  overlap, 
de-excite:  to  lose  energy  as  in  a  nucleus  changing  from  a  higher 

excited  nuclear  state  to  a  lower  one. 

Doppler  shift:  the  change  of  frequency  of  an  oscillating  wave  due 
to  motion  of  the  source,  observer  or  both, 
electric  multipole  :  A  single  point  charge  is  a  monopole.  A  dipole 
is  two  opposite  charges  separated  by  a  small  distance  s.  A 
quadrupole  is  obtained  by  displacing  a  dipole  a  distance  s^ 
and  placing  a  dipole  (separation  distance  s^)  of  opposite 
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sign  in  the  original  dipole  position  such  that  s^  =  s This 
is  continued  for  higher  poles. 

electric  multipole  radiation:  radiation  of  the  form  given  off  by  an 
oscillating  electric  multipole. 

excited  nucleus:  a  nucleus  whose  energy  is  above  its  lowest  possible 
value.  The  energies  which  can  be  assumed  by  the  nucleus  are 
discrete  up  to  the  continuum  (discrete  states) . 

gamma  ray:  highly  energetic  electromagnetic  radiation  (wavelength 
£  lO"10  m). 


Ge(Li)  detector:  lithium  drifted  germanium  solid  state  depletion  region 
detector  possessing  excellent  energy  resolution  (3~4  keV) . 

Hamiltonian:  quantum  mechanical  expression  for  the  total  energy  of  a 

system. 

magnetic  multipole:  a  concept  arising  from  the  scalar  magnetic  poten¬ 
tial 


V 


m ' 


=  Itt\ 

i  477 


where  I  is  the  current  and  $"2  the  solid  angle  subtended  by  the 
circuit  at  the  field  point  at  which  a  magnetic  induction,  _B, 
is  to  be  measured, 


Then 


m  4tt 


M*  r 


where  t  is  a  small  volume  and  Mdx  is  the  magnetic  dipole, 
magnetic  multipole  radiation:  radiation  of  the  form  given  off  by  an 
oscillating  magnetic  multipole. 
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mean  nuclear  lifetime!  Tm»  the  time  taken  for  a  number  of  excited  nuclei 

to  decay  so  that  1/e  of  the  original  number  of  excited  nuclei 

still  survive  in  the  excited  state. 

-9 

nanosecond:  nsec,  10  second. 

nuclear  half! if e:  t^ie  taken  for  a  number  of  excited  nuclei  to 

decay  so  that  1/2  of  the  original  number  of  excited  nuclei  still 

survive  in  the  excited  state.  T- =  0.693147  T 

1/2  m 

nuclear  model:  a  mathematic-physical  concept  of  the  nucleus  which 

attempts  to  present  as  accurate  a  picture  as  possible  of  the 

true  nuclear  form  (which  is  as  yet  uncertain) . 

nuclear  state:  either  the  lowest  (ground)  or  one  of  the  excited  energy 

levels  of  a  nucleus  (cf.  excited  nucleus)  characterized  by 

various  quantum  numbers  (spin,  parity,  isotopic  spin). 

parent,  daughter  nucleus:  When  a  nucleus  gives  off  a  particle  and 

becomes  a  different  nucleus,  the  original  is  termed  the  parent 

and  the  resulting  nucleus  is  termed  the  daughter. 

parity:  the  parity  describes  the  effect  on  a  system  by  transposing  each 

coordinate  through  the  space  origin.  If  the  system  does  not 

change  sign  it  is  said  to  have  even  parity;  if  it  does  change 

sign  it  is  said  to  have  odd  parity. 

photomultiplier  tube:  an  electronic  device  which  gives  a  relatively 

large  electrical  output  for  a  very  small  light  input. 

-12 

picosecond:  psec,  10  second. 

scintillator:  any  material  which  emits  light  when  struck  by  an  atomic 

or  subatomic  particle,  or  by  electromagnetic  radiation. 
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target:  the  collection  of  nuclei  to  be  bombarded  in  a  nuclear  physics 

experiment.  Usually  in  solid  or  gaseous  form. 

Time  to  Amplitude  Converter:  TAC,  an  electronic  device  which  outputs 
a  signal  proportional  to  the  time  between  the  input  start  and 
input  stop  pulses. 

timing  jitter:  uncertainty  in  timing  due  to  various  effects, 
transition  probability:  the  probability  for  a  nucleus  to  change  from 
one  excited  nuclear  state  to  another, 
walk:  change  in  the  timing  signal  due  to  change  of  the  triggering 

crossover  intercept  with  pulse  height  variation, 
wavef unction :  a  mathematical  description  of  a  given  state  of  a  quantum 


mechanical  system. 
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